Link overlaps at Criticality and Universality in Ising Spin Glasses 
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Extensive simulations are made of link and spin overlaps in four and five dimensional Ising Spin 
Glasses (ISGs). Moments and moment ratios of the mean link overlap distributions (the variance, 
the kurtosis and the skewness) show clear critical behavior around the ISG ordering temperature. 
The link overlap measurements can be used to identify the ISG transition accurately; the link overlap 
is often a more efficient tool in this context than the spin overlap because the link overlap inter- 
sample variability is much weaker. Once the transition temperature is accurately established, critical 
exponents can be readily estimated by extrapolating measurements made in the thermodynamic 
limit regime. The data show that the bimodal and Gaussian spin glass susceptibility exponents 7 
are different from each other, both in dimension 5 and in dimension 4. Hence ISG critical exponents 
are not universal in a given dimension, but depend on the form of the interaction distribution. 

PACS numbers: 75.50.Lk, 05.50.+q, 64.60.Cn, 75.40.Cx 



INTRODUCTION 

Wc have studied the equilibrium link and spin over- 
lap distributions (defined below, Eqs. ^ and ([2])) in 
some detail for Ising Spin Glasses (ISGs) on [hyperjeubie 
lattices with bimodal and Gaussian near neighbor inter- 
action distributions in dimension five and with bimodal 
near neighbor interactions in dimension four. 

The Hamiltonian is as usual 



n 



JijSiSj 



(1) 



with the near neighbor symmetric bimodal (± J) or Gaus- 
sian interaction distributions normalized to (Jfj) = 1. 
Throughout we will quote inverse temperatures /3 ~ 1/T. 

The link overlap parameter [l[ in ISG numerical sim- 
ulations is the bond analogue of the intensively studied 
spin overlap. In both cases two replicas (copies) A and 
B of the same physical system are first generated and 
equihbrated; updating is then continued and the "over- 
laps" between the two replicas are recorded over long 
time intervals. The spin overlap at any instant t corre- 
sponds to the fraction q{t) of spins in A and B having 
the same orientation (both up or both down), and the 
normalized overall distribution over time is written P{q). 
The link overlap corresponds to the fraction qe(t) of links 
(or bonds or edges) between spins which are either both 
satisfied or both dissatisfied in the two replicas; the nor- 
malized overall distribution over time is written Q{qi). 
The explicit definitions are for the spin overlap 
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and for the link overlap 



(2) 



(3) 



where N is the number of spins per sample and the 
number of links; near neighbor spins i and j are linked, 
as denoted by ij. We will indicate means taken over 
time for a given sample by (■ • ■ ) and means over sets of 
samples by [■■•]. The physical distinction between the 
information obtained from P{q) and Q{qe.) is frequently 
illustrated in terms of a low temperature domain picture 
■ " Overlap equivalence" has been proved in 0] . 

The critical behavior of the link overlaps in ISGs has 
not been studied before, as far as we are aware. It turns 
out that the moments and moment ratios of the link 
overlap distributions have characteristic forms as func- 
tions of temperature around /3c and that these data can 
be used to supplement and improve on information from 
spin overlap measurements. 

It is widely considered to be self-evident that just as 
the standard universality rules hold exactly at ferromag- 
netic ordering transitions, they should hold also for ISG 
transitions. Early exponent estimates were erratic (see 
a summary in Ref. [^). Recent careful comparisons be- 
tween the critical exponent estimates for bimodal and 
Gaussian interaction distribution ISGs have indeed con- 
cluded that both in 3d [3, and in 4d 0] the exponents 
for the two systems are the same to within numerical 
precision. High Temperature Series Expansion (HTSE) 
analyses also concluded that the estimates for the expo- 
nent 7 for different ISGs were compatible with universal- 
ity to within the precision of the method Q- However, 
the application of the standard universality rules to ISGs 
has been questioned on the basis of dynamic simulations 
[1, Q and it is relevant that the critical exponents in 3d 
Heisenberg spin glasses have been shown experimentally 
to depend on the strength of the Dzyaloshinsky-Moriya 
anisotropy [lo| . and so are not universal. 

The present simulation results concern ISGs in the high 
dimensions d = 5 and d = 4, where independent infor- 
mation from HTSE can be used in conjunction with the 
numerical data. Once accurate values of ordering tem- 
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peratures have been obtained using information from a 
combination of spin overlap, link overlap and HTSE data, 
the critical exponent 7 can readily be estimated from the 
temperature variation of the spin glass susceptibility in 
the paramagnetic state. The data presented below show 
that in the same dimension, bimodal and Gaussian ISGs 
have different values for 7, so the standard simple uni- 
versality rules are not obeyed. 



SPIN OVERLAP PARAMETERS 

The spin glass susceptibility is defined by 

x{P,L)^NMtn (4) 

where q{t) is the spin overlap Eq. ([2]). The standard 
Binder cumulant criterion which is widely used to esti- 
mate the ordering temperature /3c in ISGs consists of the 
observation of intersections of mean spin overlap kurto- 
sis curves as functions of temperature for different sam- 
ple sizes L. We use "P kurtosis" Eqn. ^ to specify the 
kurtosis of the spin overlap distribution to distinguish it 
from the kurtosis of the link overlap distribution, the " Q 
kurtosis" Eqn. A mean P kurtosis can be defined 

cither as 



Pk{fi,L) 



id*) 



which we will use here, or alternatively 



Pkm{f3,L) = 



Ml 



(5) 
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which is the definition more frequently used. P kurtosis 
data are generally expressed in terms of the Binder cu- 
mulant g{P,L) = [3 — Pkm{f3, L)]/2. One drawback to 
this procedure for estimating f3c is that the inter-sample 
variability of (g^) and a fortiori of (g'*) are strong at f3c in 
ISGs. The normalized variation of the ISG susceptibility 
(the non-self-avcraging parameter) 
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(7) 



is typically about 0.20 at /3c [1, [TH [l^, and results on 
large numbers of samples must be recorded at each size to 
overcome statistical fluctuations in the mean P kurtosis. 
There are in general finite size corrections so 



Pfe(/3c,L) = Pfe(/3c,oo)(l + CfcL- 



(8) 



with prefactor Ck and exponent u which are a priori un- 
known. Delicate extrapolation to infinite L is required to 
estimate the thermodynamic limit critical temperature. 
To obtain the intersection point of curves Pk{/3,L) and 
Pk{/3, L ) it is necessary to equilibrate samples up to the 
larger size L while the position of the intersection point 



is affected mainly by the larger finite size correction of 
the smaller size L. 

As well as higher order moment ratios, there are other 
dimensionless parameters of the spin overlap distribu- 
tions which also can be studied, if the one-sided distri- 
butions of the absolute value of the spin overlap Pabsd^l) 
are recorded. These include the second moment ratio 



Pi 



w 



(9) 



the mean skewness of the absolute spin overlap distribu- 
tion 



absskcw 



((kl-(kl))-^) 
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and the mean kurtosis of the absolute spin overlap dis- 
tribution 



P 



abskurt 



(11) 



It turns out that the first two moment ratio parame- 
ters are of the standard phenomenological coupling form; 
at /3 = each takes up a standard value corresponding to 
that of a one-sided Gaussian, and then decreases towards 
a low value at high (3. As functions of /? at fixed L, at /3c 
to leading order the curves go through a size- independent 
critical value and have a maximum slope [dP(/3, L)/d/3]c 
whose value increases as L'^^'^ . The parameters generally 
show finite size corrections. On the other hand data on 
the novel parameter Pabskurt(/3, i) defined by Eqn. (fTTj) 
shows a deep dip as a function of temperature for each L. 
With increasing L the dip narrows and the position of the 
minimum is tending towards /3c subject to a weak finite 
size correction. Pabskurt(/3, i) behaves in quite different 
ways in the ISGs and in the ferromagnet [13[. This spin 
overlap based parameter is obviously a useful supplemen- 
tary measurement for estimating critical temperatures. 

We will not discuss here the correlation length ratio 
^{p,L)/L which is a widely used phenomenological cou- 
pling parameter independent of the spin overlap distri- 
bution. 



LINK OVERLAPS 

Turning to the link overlaps, for the Gaussian ISG it 
has been shown IJ] that in equilibrium 



P(/3) 



(l-[(|C/|)(/3,L)]//3 



1 



(12) 



where [(|t/|)(/3, L)] is the mean energy per bond. The 
5d Gaussian ISG data presented here satisfy this equi- 
librium condition over the full temperature range used; 
the bimodal samples equilibrated faster than the Gaus- 
sian ones and were equilibrated for as long times so we 
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will consider that for present purposes effective equilibra- 
tion has usually been reached. However, the data show 
that the condition Eq. ([T^ is necessary but is not strin- 
gent enough to guarantee true equilibration. A stricter 
and more general condition, which can be applied what- 
ever the interaction distribution, is that all spin over- 
lap or link overlap parameters for each individual sample 
should vary smoothly with temperature. By inspection 
of the individual sample data sets it can be seen if and 
when the equilibrium condition begins to break down as 
the temperature is lowered. This test shows that for the 
same L some sam ples equilibrate more easily than others, 
as noted in Ref. [l5|. What has not been remarked on 
is that "simple" samples, where the spin overlap distri- 
bution is tending to two pure peaks beyond the ordering 
transition and (qi) is high, equilibrate more easily than 
"complex" samples, those for which the spin overlap dis- 
tribution remains multi-peaked and (qi) is low even at 
low temperatures. 

For the symmetric bimodal ISG there are simple rules 
on the mean link overlap (qi). If Ps{P) is the probability 
that a bond is satisfied, by definition 



|C/(/3)| = l-2p,(/3) 



(13) 



and for symmetric interaction distributions where the 
Nishimori point is at /3 = (uncorrelated satisfied bond 
positions) a strict lower limit on {qi) is given by 

[{qi)iP,L)]>pl + (l-ps)^ - 2ps{l-ps) = [[/(/?,£)"]. 

(14) 

In the high temperature limit \U\{I3) tanh(/3) so 
i?(/3) — 3. As /3 increases R{P) drops and appears to 
tend gradually towards 1. 

For a pure near neighbor ferromagnet (so with transla- 
tional invariance) [(g^ )(/3, L)]/[J7(/3, L)^] = 1 at all tem- 
peratures [l^. For the bimodal ISG this ratio is equal 
to 1 for small /3 but then gradually grows as f3 increases 
and certain bonds become preferentially satisfied. 

Certain moments and moment ratios for the link over- 
lap distributions Qi^qi) show characteristic critical be- 
havior in the ISGs, as they do in a ferromagnet [l3l |. 
Link overlap data for a given (3 and L can be recorded 
for virtually no extra computational cost in a simula- 
tion designed for spin overlap measurements, while the 
inter-sample variability of the link overlap distributions 
at Pc is considerably weaker than that of the spin over- 
lap distributions. This implies that measurements on 
mean link overlap values require far fewer samples than 
those on mean spin overlap values, or alternatively that 
with the same number of samples the mean link over- 
lap measurements are more precise than the mean spin 
overlap measurements. Thus link overlap based data are 
more efficient for obtaining accurate estimates of critical 
temperatures, and so of critical exponents, than are spin 
overlap data. 



As well as the mean link overlap [{qi)], and the variance 
of the link overlap 



(15) 



we have recorded three dimensionless moment ratios for 
each sample and their averages over all samples. These 
are the Q kurtosis 



QkW,L) 

the Q skewness 



{m - {<ii)t 
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iqi~{qi)y 
(qe - {qe)y 



3/2 



(17) 



and Qw which is the mean squared signal to noise ratio, 
or the mean of the inverse square of the coefficient of vari- 
ation. This has a clumsy name but a simple definition: 



{qi)' 



Nidqe ~ {qem 



(18) 



As we will see, its derivative with respect to /3 has a 
minimum which location is very close to /3c. Finally, the 
quantity Qv is defined as the squared ratio of the mean 
deviation and the standard deviation (thus a relative of 
Pw above) i.e. 



{\qi-{qi)\)' 



{{qi (qe))') 



(19) 



Obviously can only be found by storing the actual dis- 
tributions Q{qi) for individual samples during simulation 
rather than just the raw moments. 

At high temperatures, i.e. /3 — )> 0, for both bimodal 
and Gaussian ISG interactions the Q{qi) distributions 
become symmetric, Gaussian, and centered on qi = 0, 
so Qfc(O) = 3,Qs(0) = 0,Q,„(0) = 0. As /3 is increased 
through /Sc the Q{qi) distributions become fat tailed and 
asymmetric, so Qk{P, L) and Qs{l3, L) show peaks in the 
region of Pc- 

The amplitudes of the peaks decrease with increas- 
ing L; this could be called an "evanescent" critical phe- 
nomenon as it will disappear in the thermodynamic limit. 
Allowing for a weak finite size correction term, the po- 
sitions of the maxima /3max for each set of peaks tend 
towards /3c with increasing L. 

Physically, the peaks in the excess Q kurtosis ("fat 
tailed" distributions) and the Q skewness near /3c in ISGs 
must be related to the build up of inhomogeneous tempo- 
rary correlated spin clusters around criticality. The data 
show that these clusters do not produce a visible effect 
on the form of the Q{qi) distribution until L is smaller 
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than the thermodynamic correlation length £,{(3), which 
is related to the typical cluster size. For larger L the 
cluster effects average out in Q{qe). Only when the ratio 
L/£^{/3) is smaller than some value so that a cluster can 
englobe the entire sample do deviations from the Gaus- 
sian form appear. As ^(/3) diverges at /3c, the Q kurtosis 
and the Q skewness will each tend to a peak for fixed L, 
and the peaks can be expected to be situated exactly at 
/3c in the large L limit. Indeed analogous behavior can 
be observed in a pure Ising ferromagnet, with excess Q 
kurtosis and Q skewness peak positions tending towards 
/3c with increasing L [l^. 

The Qvar(/3,i) and Qw{l3,L) parameters are closely 
related as [{qi){(3,L)] becomes almost independent of L 
at large L. We will show data for log((5var(/3, L) — 1) and 
Qw(J3,L). Both of these parameters have basically the 
form of a phenomenological coupling, with curves for dif- 
ferent L intersecting at crossing points /3cross(^, L ) which 
approach /3c as L and L are increased. The corrections 
to scaling for the two parameters are slightly different. 
As the inter-sample variability for these parameters is 
much weaker in the region of /3c than is the equivalent 
variability for the phenomenological couplings based on 
spin overlap distributions, they provide an accurate tool 
for estimations of the critical temperature. 



CRITICAL EXPONENT ESTIMATES 

The thermodynamic limit spin glass susceptibility in- 
cluding the Wegner confluent correction to scaling term 
M is 



(20) 



with the natural ISG scaling variable t = 1 — (/3//3c)^ 
l[i3|. The analo gous correlation length expression is 



^(r)//3 = C^T-'il + ay +Kt + ■■■). 



(21) 



with the same exponent 9. (Unfortunately ISG suscep- 
tibility and correlation length data are most frequently 
analyzed using the scaling variable t = {T — Tc)/Tc = 
/3c//3 — 1 which is inappropriate for ISGs except as an ap- 
proximation in a narrow region around the critical point.) 
Once the value of /3c can be considered to be precisely de- 
termined, a further step is to make a plot of the temper- 
ature dependent effective exponent j{t) from the 
spin glass susceptibility data x(/3), with the definition 



7(r) = -dlog(x(/?))/dlog(r) 



(22) 



and the temperature dependent effective exponent i^(t) 
can be defined by 



In a [hyperjcubic lattice, from the second term in the ISG 
HTSE the limiting effective exponent at infinite temper- 
ature 7(1) = 2(i/3^ exactly, where 2d is the number of 
nearest neighbors. The ratio 7(0)/7(l) is directly related 
to the strength and sign of the confluent correction coefH- 
cient . If the remaining correction terms are negligible, 
then to leading order 



7(0) - 2d/32 



(24) 



which gives a criterion for the strength and sign of the 
confluent correction to scaling term (but not for the value 
of the exponent 6) once /3c and 7(0) have been estimated. 

As long as L 3> ^(/3) the finite size numerical data are L 
independent and so are effectively in the thermodynamic 
limit infinite size regime. The region where this condition 
holds for each L can be seen by inspection of 7eff (t, L) 
and other parameters such as foff(T, L). Extrapolating 
the thermodynamic limit regime susceptibility data to 
criticality at r = provides an estimate for the critical 
7(0), see e.g. Fig. [181 This value is fully reliable under 
the unique condition that /3c has been correctly deter- 
mined. An accurate knowledge of /?c is essential; there 
is a one-to-one relationship between the estimate for the 
critical 7(0) and the value of /3c taken to construct the 
plot. Finally, with fixed /3c and 7(0) the scaled spin glass 
susceptibility can be plotted in the form xi''')''''' against 
where the correction to scaling exponent 9 is chosen 
such that the plot is a straight line over the thermody- 
namic limit data region. The parameters and a-^ can 
be read off this plot. 

It can be noted that in the thermodynamic limit regime 
L > ^(r) there is "self-averaging", or in other words all 
individual ISG samples of a system have the same prop- 
erties and in particular the same spin glass susceptibility. 
Thus in this regime there is no real need to average over 
large numbers of samples to obtain accurate measure- 
ments of the mean x(t). In addition, equilibration in the 
thermodynamic limit region is relatively rapid so mea- 
surements are very reliable and not subject to equilibra- 
tion difficulties. On the contrary, in the regime near, at 
and beyond /3c "lack of self-averaging" sets in; the inter- 
sample variability is important. The non-self-averaging 



i^ir) = -dlog(^(/3)//3)/dlog(r) 



(23) 



parameters are size independent jllL Il2| at and beyond 
/3c so however large the individual ISG samples they are 
all different from each other; there is a wide distribution 
of values of the spin glass susceptibility and other prop- 
erties. The onset of a non-zero variability is a spin glass 
criterion for an approach to the transition temperature 
which obviously has no equivalent in pure systems such 
as simple ferromagnets. Even in diluted ferromagnets 
the non-sclf-averaging is non-zero in the thermodynamic 
limit only at /3 = /3c. 
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NUMERICAL SIMULATIONS 

For equilibration and measurements we used standard 
heat bath updating (without parallel tempering) on ran- 
domly selected sites. The samples (usually 64) started 
off at infinite temperature and was then gradually cooled 
before reaching their final designated temperature. For 
temperatures near Tc this means that each sample went 
through at least 10^, sometimes 10^, sweeps before any 
measurements took place. Normally there were about 
10 sweeps between measurements, depending on tem- 
perature, maintaining on average X'* spin flips between 
each measurement. For each sample and temperature we 
collected between 10^ and 10^ measurements depending 
on lattice size. The test for equilibration was discussed 
above. It can be noted that a sample with L = 8 in di- 
mension 5 corresponds to as many individual spins as a 
sample with L = 32 in dimension 3. 

BIMODAL ISG IN DIMENSION FIVE 

No detailed numerical simulation measurements have 
been reported before for ISGs in dimension five. How- 
ever, the analysis of a 15 term High Temperature Se- 
ries Expansion (HTSE) calculation [7[ gave the estimate 
Pi = 0.154(3), i.e. (3^ = 0.3925(35), with a critical expo- 
nent 7 = 1.95(15), for the bimodal ISG in 5d. We have 
re-analyzed the two series in an unorthodox but transpar- 
ent manner, Appendix I, and obtain values for (3^ and 7 
very close to the central values in the original HTSE anal- 
ysis but with additional information as to the strength of 
the correction to scaling term in the two cases. It can be 
seen that the correction to scaling is strong in the 5d bi- 
modal case (and practically negligible in the 5d Gaussian 
case). 

Numerical data derived from the various Pq{(3,L) and 
Qgj, (/3, L) distributions, all taken in the same runs on the 
same sets of samples for each L, are shown in Figures[l]to 
[71 The error bars correspond to inter-sample variability 
for each particular parameter. 

The spin overlap P{q) based phenomcnological cou- 
plings P kurtosis Eq. ([5]), Pw Eq. and the skewness 
of the absolute P distribution Eq. (|10p show very simi- 
lar forms. Because of the strong dispersion of individual 
sample parameter values, crossing points derived from 
the present data with a modest number of samples scat- 
ter and cannot provide an accurate estimate for /3c. The 
data are broadly compatible with the central HTSE es- 
timate but on their own do not provide anything like a 
critical test, which for these parameters would require 
averaging over a much larger number of samples. 

With the present sets of of samples, the most use- 
ful spin overlap parameter is the kurtosis of the abso- 
lute value distribution, Pabskurt(/3, i) Eq. (fTTI) . It shows 
a strong dip for each L in the region of the HTSE Pc 



estimate. The center of the dip as a function of tem- 
perature can be estimated quite accurately for each L. 
When the central dip positions are plotted against 1/L 
and extrapolated to 1/L = 0, the data give an estimate 
Pc = 0.392(3). This is the most precise estimate obtained 
from the spin overlap distributions and is fully consistent 
with the HTSE central value. 

The link overlap parameters Qyj{f3,L) and log((5var — 
l)(/3, L) show much smaller inter-sample variability than 
the parameters based on spin overlap. By luck log(Qvar — 
l)(/3, L) has negligible finite size correction to scaling for 
the 5d bimodal case; the curves for different L all inter- 
sect at the same /3cross(-^, L) = 0.392(1), which coincides 
with the central value for (3c from the HTSE estimate 
13c = 0.3925(35) 0. This result both vahdates the as- 
sumption that this link overlap parameter is a bona fide 
phenomcnological coupling, and improves the precision 
on the value of the critical temperature. For Qy^{f3,L) 
there are weak finite size corrections to scaling, but the 
intersection points /3cross{L,L) as functions of L for 
fixed L can be extrapolated to infinite L to obtain an 
estimate /3c = 0.392(3) which is again in agreement with 
the HTSE value. In addition the (negative) maximum in 
the derivative dQ^iP, L)/df3 deepens with increasing L 
and its position f3nia,x{L) can with extrapolation also be 
used to estimate (3c- The intersection criterion and the 
maximum slope criterion conveniently bracket the criti- 
cal Pc more and more closely as the sizes are increased. 
From these Qvar and Qu, numerical data alone one can 
thus derive a very precise estimate Pc = 0.392(1), entirely 
consistent with the HTSE central estimate. The Q kur- 
tosis Qk and Q skewness Qs peak positions are subject 
to finite size corrections but are also consistent with this 
estimate for f3c- 

Thus the 5d bimodal ISG is a particularly favorable 
case to validate the assumption that the link overlap pa- 
rameters show strictly critical forms in ISGs as they do 
in a ferromagnet [ist . 

THE GAUSSIAN ISG IN DIMENSION 5 

The critical temperature was estimated from the HTSE 
analysis to correspond to P^ = 0.174, 0.176(3)or 0.177(3) 
according to the different analysis techniques [7|, i.e. Pc = 
0.420(3). The spin glass critical exponent was estimated 
by the HTSE analysis to be 7 = 1.75(15). 

The data for spin and link overlap moments and mo- 
ment ratios are shown in Figs [8] to [TlJ The general form 
for each parameter is similar to that for the 5d bimodal 
ISG, with the appropriate Pc = 0.421(2) being estimated 
from the absolute P kurtosis and from Q^, allowing for 
corrections to scaling, in agreement with the central value 
from the HTSE analysis. 

In the data for for the Q kurtosis and the Q skew- 
ness. Fig. [TU] and Fig. [TTl the peak positions are tending 
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FIG. 1: (Color online) The mean P kurtosis Eq. [5] for 5d 
bimodal interaction samples with L = 4 (blue circles), L — 
5 (green squares), L — 6 (red triangles) and L = 8 (black 
diamonds). The vertical red line corresponds to the HTSE /3c 
central value 



FIG. 3: (Color online) The mean kurtosis of the absolute value 
spin overlap distribution, Pabskurt(/3, i), Eg. Illl for the sets 
of 5d bimodal interaction samples L = 4, 5, 6, 8 (sizes coded 
as in Fig. [T]). The vertical red line corresponds to the HTSE 
/3c central value. 




FIG. 2: (Color online) The mean Pw parameter Eq. [9] for 
the sets of 5d bimodal interaction samples with L = 4, 5, 6, 8 
(sizes coded as in Fig. [T]). The vertical red line corresponds 
to the HTSE /3c central value. 

towards /3c with increasing L more slowly than in the bi- 
modal case. This may be due to strong finite size correc- 
tions or possibly a peculiarity of the Gaussian interaction 
distribution. 

THE BIMODAL ISG IN DIMENSION 4 

From an analysis of HTSE data for the 4d bimodal 
ISG, Daboul et al estimate /S^ = 0.26(2), i.e. Pc = 
0.51(2). (HTSE estimates in 4d are intrinsically less pre- 
cise than in 5d). A critical temperature /3c = 0.493(7) 
was estimated [l9| from simulation measurements of high 



FIG. 4: (Color online) The logarithm of the Q variance Eg. 1151 
minus 1, log((3var — 1)(/3, L), for the sets of 5d bimodal interac- 
tion samples with L = 4, 5, 6, 8 (sizes coded as in Fig.[T]). The 
vertical red line corresponds to the HTSE /3c central value. 
The horizontal red line represents the HTSE /?c range 0]. 

statistical accuracy to L 10 using the Binder parameter 
crossing point criterion, but corrections to scaling were 
not allowed for. A further estimate is /3c = 0.5025(25) 
[2^ from unpublished Binder parameter data to L = 12 
by A. P. Young. From extensive domain wall free en- 
ergy measurements to L = 10 Hukushima gives an es- 
timate Tc = 2.00(4) [iOl, i.e. /?c = 0.50(1). In fact 
the raw data show significant finite size corrections, 
which affect the extrapolated estimate for the value of 
/3c in the infinite size limit. This can be seen clearly 
in the data shown in Fig. 4 of Ref. the cross- 

ing points Tcross{L, L*) evolve regularly from the small- 
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FIG. 5: (Color online) The inverse normalized Q variance 
Qw{l3, L) Eq. 1181 for the sets of 5d bimodal interaction sam- 
ples with L = 4, 5, 6, 8 (sizes coded as in Fig. [T]). The vertical 
red line corresponds to the HTSE /3c central value. 



FIG. 7: (Color online) The Q skewness Eqn. Qs{l3,L), 
for the sets of 5d bimodal interaction samples with L — 
4,5,6,8 (sizes coded as in Fig. The vertical red line cor- 
responds to the HTSE /3c central value. 
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FIG. 6: (Color online) The Q kurtosis Eq. [161 Qk{l3,L), for 
the sets of 5d bimodal interaction samples with L = 4, 5, 6, 8 
(sizes coded as in Fig. [T|). The vertical red line corresponds 
to the HTSE /3c central value. 



est sizes to the largest sizes measured : Tcross(4, 6) 



2.20,r„oss(6,8) 



2.06,Tcross(8,10) ~ 2.00. By inspec- 
tion, the infinite size limit crossing temperature must be 
distinctly lower than T = 2.00, i.e. /3c > 0.50 [U. 

Simulations were carried out on sets of 64 samples of 
size L = 3, 4, 5, 6, 7, 8, 9 and 12. As in 5d the phenomeno- 
logical couplings based on the spin overlap were strongly 
affected by the inter-sample variability; much larger sets 
would have been needed to obtain crossing point data 
for these parameters of similar statistical precision as in 
the earlier results for the Binder parameter. However, 
the position of /3dip, the minimum of the dip in Pabskurt, 
FigHH is independent of L at /3dip = 0.506(3) to within 
the statistical precision. 



FIG. 8: (Color online) The Pabskurt(/3, L), Eq.[TTl for the sets 
of 5d Gaussian interaction samples with L — 4,5, 6, 7, 8, 10 
(blue circles, green inverted triangles, red triangles, olive 
squares, black diamonds, pink left triangles) . The vertical 
red line corresponds to the HTSE central value /3c = 0.421. 



For the link overlap parameter Qy^{f3,L), Fig. [T51 and 
Fig. I14[ the inter-sample variability is much weaker than 
for the spin overlap phenomenological coupling parame- 
ters, so the crossing points for successive L are better de- 
termined, but there are both finite size corrections with 
the crossing points evolving towards larger f3cross{L, L* 
with increasing L,L*, and odd-even effects in L. The 
/3cross(i — 1,-^ + 1) values for crossing points between 
Q.u]{(3,L — 1) and Qu,{f3,L + 1) are shown against 1/L 
in Fig. [13 Extrapolating to 1/L = 0, /3c = 0.510(5) 
from this criterion. The deviation ratio Qy, Fig. 1161 also 
show weak sample variability with a distinct minimum 



FIG. 9: (Color online) The Qn,(P,L) parameter of Eq. [HI 
for the sets of 5d Gaussian interaction samples with L = 
4,5,6,7,8 (sizes coded as in Fig. |8]). The vertical red line 
corresponds to the HTSE central value /3c = 0.421. 
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FIG. 10: (Color online) The Q kurtosis Qk{l3,L) Eq. [H 
for the sets of 5d Gaussian interaction samples with L = 
3,4,5,6,7,8,10 (sizes coded as in Fig. |8]). The vertical red 
line corresponds to the HTSE central value/3c = 0.421. 

approaching /3c. The Q kurtosis Qk{(i,L) increasing L 
towards a limiting /? values for 1/X = which is con- 
sistent with the estimate from Q^j, Fig. [T71 It can be 
concluded that /3c = 0.510(5) in full agreement with the 
central HTSE estimate and with the previous numer- 
ical measurements once finite size corrections are fully 
allowed for. 



THE CRITICAL EXPONENT 7 

Daboul et al concluded that in each dimension the 
HTSE critical bimodal and Gaussian 7 values for the dif- 
ferent interaction distributions which they studied were 
compatible with universality in ISGs to within the un- 



P 



FIG. 11: (Color online) The Q skewness Qsil3,L) Eg. [TTl 
for the sets of 5d Gaussian interaction samples with L — 
3, 4, 5, 6, 7, 8, 10 (cyan stars, blue circles, green inverted tri- 
angles, red triangles, olive squares, black diamonds, pink left 
triangles) . The vertical red line corresponds to the HTSE 
central value /3c = 0.421. 
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FIG. 12: (Color online) The absolute P kurtosis 
fabskurt(/3, i), Equ. (|lip . for the sets of 4d bimodal inter- 
action samples with L = 3, 4, 5, 6, 7, 8, 9, 12 (cyan stars, blue 
circles, green inverted triangles, red triangles, olive squares, 
black diamonds, purple right triangles, pink left triangles). 
The vertical red line corresponds to the HTSE central value 
/3c = 0.51 



certainties of the HTSE analysis. However, their error 
bars for each 7 value were relatively large, as they did 
not have access to simulation data which supplement the 
HTSE calculations and which refine both the /3c and the 
7 estimates. 

The effective 7(r, L) values are defined by Eq. (|22p . 
see the "Critical exponent estimates" section. For the 
5d bimodal and 5d Gaussian systems 7(r, L) is shown 
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FIG. 13: (Color online) The parameter Q^{I3,L), Eqn. dHJ, 
for odd L sets of 4d bimodal interaction samples with L — 
3, 5, 7, 9 (size coding as in Fig. I12|) . The vertical red line cor- 
responds to the HTSE central value /3c ~ 0.51 



FIG. 15: (Color online) The Q kurtosis Qk{P,L), Eqn. JTS}, 
for sets of 4d bimodal interaction samples with L = 4, 5, 6, 7, 8 
and 12 (size coding as in Fig. I12|l . The vertical red line cor- 
responds to the HTSE central value /3c = 0.51 




FIG. 14: (Color online) The parameter Q„(/3,L), Eqn. ([1811 . 
for even L sets of 4d bimodal interaction samples with L = 
4,6,8,12 (size coding as in Fig. I12|) . The vertical red line 
corresponds to the HTSE central value /3c = 0.51 



in Fig. [18] and Fig. [TH with in each case r = 1 — 
[PI PcY, the data being plotted using the optimal val- 
ues /3c(bimodal) = 0.392 and /3c(Gaussian) ~ 0.421 as 
quoted above. See also Figs. [20l and [211 for a plots of the 
reduced susceptibility in these cases with the assumed 
values of /3c and critical exponents. 

It can be noted that in 5d a sample with L = 8 is 
"large" in the sense that for this size xi'TjL) remains in 
the thermodynamic limit until r 0.2. To approach 
/3//3c = 1 equally closely in 3d would require samples 
with L ~ 100. The difference is a consequence of the 
much higher value of the exponent v in dimension 3, 
which is not far from the lower critical dimension where 
v diverges. 



FIG. 16: (Color online) The Q deviation ratio Qv{l3,L), 
Eqn. (|19|l . for sets of 4d bimodal interaction samples with 
L = 3, 4, 5, 6, 7, 8, 9 and 12 (size coding as in Fig. [12)). The ver- 
tical red line corresponds to the HTSE central value /3c = 0.51. 

For the 5d bimodal ISG the plot of the effective j{t) 
shows a critical limit for the extrapolated thermodynamic 
limit regime 7(0) = 1.92(5), and a strong correction to 
scaling with exponent 9^1. With this value of 7 in hand 
x(t)t^ can be plotted as a function of r, Fig. [19l which 
can be fitted by = 0.59(1 + 0.69t). So finally 

for the 5d bimodal ISG, in Eq. ([101) 7 = 1.92(5), 6* - 
1,C^ = 0.59, ajj- = 0.69. As ~ 1 the second correction 
term in Eq. ([^(I)) cannot be distinguished from the leading 
term. The fit curve for 7eff(T) in Fig.[T9|is the derivative 
of xi''') with the same parameters. The parameter set 
including the sign and the approximate strength of 
are consistent with those which can be estimated entirely 
independently from the HTSE data (see Appendix I). 
The critical 7 and values are consistent with but more 
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FIG. 17: (Color online) The Q kurtosis Qk{P,L) peak tem- 
perature /3niax(i) (black squares), the Q-u minimum location 
(blue circles) and the AQw/Aji minimum location (red trian- 
gles) against 1/ for sets of 4d bimodal interaction samples 
with L = 3, 4, 5, 6, 7, 8, 9. See Figs. [H [H [15] and [16] (deriva- 
tive). The extrapolation to infinite L gives /3c = 0.510(5). 
Error bars were obtained from statistical resampling of half 
the data set. 



accurate than the HTSE value 7 — 1.95(15) and 6 ^ 1.1 
of Daboul et aZ ■ 

For the 5d Gaussian ISG the plot of the effective 7(t) 
shows a critical limit 7(0) = 1.66(3), = 1.03. The 
Daboul et al HTSE estimate is 7 = 1.75(15). The cor- 
rection to scaling is very weak except for a high order 
term in the region far from criticality, so no estimate 
can be made for 6. Again the HTSE data, as analysed 
in Appendix I, are in complete agreement with all these 
conclusions : a very similar critical temperature, a very 
similar critical exponent, and the same weak correction 
to critical scaling. 

The principal conclusion which can be drawn from the 
5d results, both from the simulations and from the HTSE 
data analyses, is that the critical exponent 7 — 1.92(5) 
for the bimodal ISG is significantly higher than the crit- 
ical 7 = 1.66(3) for the Gaussian ISG. 

The 4d bimodal 7(r) plot assuming /3c = 0.51 are 
shown in Fig. 1221 Extrapolating the thermodynamic 
limit regime curve to r = gives a critical exponent es- 
timate 7(0) = 3.25(10) where the error bar corresponds 
principally to the residual uncertainty in j3c. The fit pa- 
rameters to x(r)r''', Fig.[23l are = 0.30, = 2.3 and 
9 ^ 1.6, so there is a very strong correction to scaling. 
The HTSE analysis led to the estimates 7 = 2.5(3) 
and ~ 1.5. A re-analysis of unpublished 4d bimodal 
x(/3) and ^(/3) data of Ref. [2l[, are in full agreement with 
the present data as far as x(/3) is concerned. Fixing 0r. = 
0.51 and defining i/eff(r) = dlog(^(r)/^)/dlog(r) [l^, 
leads to an estimate for the correlation length critical 
exponent i' — 1.30(5). 



High precision simulation measurements have been 
made of the 4d Gaussian ISG and of a 4d diluted bimodal 
ISG 0. The critical temperature for the 4d Gaussian 
ISG was estimated from Binder parameter and correla- 
tion length ratio measurements to be /3c = 0.554(3) in 
full agreement with earlier simulation estimates 0.555(3) 



/3, 



2J1 and with the HTSE estimate P'^ = 0.314(4) i.e. 
0.560(3). For the Gaussian and the diluted bimodal 
ISG the critical exponents were estimated to be 77 = 
-0.275(25) and ly = 1.02(2) so 7 = (2 - 77)1/ = 2.32(8), 
and T] = -0.275(25), = 1.025(15) so 7 = 2.33(6) respec- 
tively. The simulation data showed that finite size correc- 
tions to scaling are weak. This is consistent with the cri- 
terion for given above. With the Gaussian critical pa- 
rameters : ~ 7-2:/3c = 2.32(8)-2.45(2) = -0.13(10); 
the Wegner susceptibility correction to scaling amplitude 
will be small. 

Naturally it was concluded in Ref. j(| that as diluted 
bimodal and Gaussian 4d ISGs have the same critical 
exponents to within the statistical uncertainties, univer- 
sality is confirmed. However, the critical 7 ~ 3.25(10) 
estimated above for the undiluted 4d bimodal ISG is sig- 
nificantly higher than the 7 ^ 2.33 values estimated for 
the Gaussian and the diluted bimodal ISGs. It happens 
that at the particular diluted bimodal bond concentra- 
tion studied in Ref. [l], p = 0.35, the kurtosis of the 
bond distribution is 1/0.35, so almost exactly equal to 
the Gaussian distribution kurtosis which is 3. It is tempt- 
ing to speculate that there could be a universality rule 
for ISGs such that at fixed dimension, the exponents de- 
pend on the kurtosis of the interaction distribution, just 
as the critical temperature of an ISG in each dimension is 
a function of the kurtosis of the interaction distribution 
25| . For the moment we have not studied the ISGs in 
dimension 3. 



CONCLUSION 

The moments and moment ratios of the link overlap 
distributions in ISGs show well defined critical proper- 
ties, analogous to those observed for the link overlap dis- 
tribution moments in a simple ferromagnet [isj . The 
inter-sample variability of the link overlap parameters in 
the ISGs is weaker than that of the spin overlap param- 
eters, so link overlap critical measurements, even with 
modest numbers of independent samples, are intrinsically 
more precise than spin overlap measurements. 

Link overlap critical data have been used here to sup- 
plement spin overlap data and HTSE analyses in order to 
obtain accurate estimates for the ordering temperatures 
of ISGs in dimensions 4 and 5. We have also introduced a 
useful spin overlap dimensionless parameter, the absolute 
P distribution kurtosis Pabskurt, Eqn. (jlip. which has not 
been previously studied. The /3c values estimated from 
these simulations are all in excellent agreement with the 
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FIG. 18: (Color online) The efltective exponent 7(r) Ea.[22lfor 
the bimodal 5d ISG samples with L = 4, 5, 6, 8 (size coding 
as in Fig. I12p . assuming /3c ~ 0.3925. The extrapolation to 
r = gives 7 = 1.92(5). 
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FIG. 19: (Color online) The effective exponent 7(r) Eq. [22] 
for Gaussian 5d ISG samples with L = 4,6, 8, 10 (size coding 
as in Fig. I12|) . assuming /3c ~ 0.421. The extrapolation to 
r = gives 7 = 1.66(3). 

entirely independent central estimates from HTSE anal- 
yses [7|, but the estimates for the critical exponent 7 
are more accurate. It should again be underlined that 
to obtain precise estimates of critical exponents it is es- 
sential to first establish reliable values for the critical 
temperatures. Once the ordering temperatures in hand, 
the effective critical exponents 7(t) can be readily and 
reliably estimated from the appropriate derivative of the 
spin glass susceptibility simulation data, Eq. (1^^ . which 
can be extrapolated to r = to obtain 7. 

The present critical 7 estimates can be compared with 
values obtained from analyses of the HTSE coefficients 
0] for both 5d and 4d, and compared with published 
simulation results Idl for the 4d Gaussian case. The sim- 




FIG. 20: (Color online) The reduced susceptibility xi''')''''' 
for the bimodal 5d ISG samples L = 4, 6, 8 (size coding as 
in Fig. assuming 13c = 0.3925, 7 = 1.92. The overall 
thermodynamic limit fit is x(t) = 0.59r~^-^^[l + 0.69r]. 
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FIG. 21: (Color online) The reduced susceptibility xi''')''''^ 
for the Gaussian 5d ISG samples L = 4, 6, 8, 10 (size coding 
as in Fig. [8|, assuming /3c = 0.421, 7 = 1.66. The overall 
thermodynamic limit fit is xi''') = 1.03r^^ '''' [1 -f • • ■]. 



ulation data show that both in 5d and 4d the critical 7 
value for the bimodal ISG is significantly higher than the 
critical 7 values for the Gaussian ISG. 

The well established universality rules which apply to 
standard second order transitions are that systems hay- 
ing the same spatial and spin dimensionalities all have 
identical critical exponents. From the present data it 
can be concluded empirically that different, more com- 
plicated, rules govern universality classes in Ising spin 
glasses. It should be remembered that ISG transitions 
are qualitatively very different from standard second or- 
der transitions. For an Ising ferromagnet in the regime 
below the Curie temperature there are just two mirror 
image families of spin up and spin down states. For an 
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FIG. 22: (Color online) The effective exponent 7(r) Ea.[22lfor 
the bimodal 4d ISG samples with L = 4,6, 8, 12 (size coding 
as in Fig. I12p . assuming /3c = 0.51. The extrapolation to 
r = gives 7 = 3.2(1) 
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APPENDIX I 

HTSE calculations in ISGs Q produce a set of terms 
for the spin glass susceptibility of the form 

Xil3^) ^ 1 + ail3^ + a2^^ + ■ ■ ■ (25) 

Each coefficient a,; is exact, but the series is in practice 
limited in length. The published ISG calculations have 
15 terms. 

This series can be compared to the mathematical iden- 
tity 




FIG. 23; (Color online) The reduced susceptibility xi''')''''' for 
the bimodal 4d ISG samples L = 4,6, 8, 12 (size coding as in 
Fig. [HI, assuming /3c = 0.51, 7 = 3.2, 6 = 1.6. The overaU 
thermodynamic limit fit is x(t) = 0.30r"3-2[l -I- 2.3r^-^]. 



ISG the non-self-averaging behavior means that at and 
beyond /3c even in the thermodynamic limit each individ- 
ual sample has different properties, in particular a differ- 
ent spin overlap distribution and so a different spin glass 
susceptibility. It is not obvious that the powerful renor- 
malization arguments which are so effective in standard 
transitions can be applied in the same manner in this 
context. It would be of interest to explore from funda- 
mental principles which relevant parameters determine 
critical exponents in the spin glass family of transitions. 



[1 — X) ' = 1 + 7X H X 



7(l + 7)(2 + 7)^3^... (26) 



for which the ratio of successive coefficients a,, is 



an 



= 1 



7-1 
n+1' 



(27) 



In the simplest case of a physical system with a critical 
spin glass susceptibility xW^) = C'x(l ~ / l^c)~'' ^ 
ratio of successive terms a„+i/a„ in Eq. (|25p would be 
(1//3^)(1 -I- (7 - l)/('^ + !))• This suggests a graphical 
analysis in terms of a plot of this ratio against l/{n + 
1), which is indeed a traditional technique for analyzing 
HTSE coefficients [3, [2^. There are two complications. 
One is the Wegner confluent correction to scaling jTij : 
the critical susceptibility is 

x(r) = C,(r)-^(l + axr« + ..-) (28) 

where r = 1 — /3^//3^. The ratios become (see jl8| ) 



1 



«x^r(7) 



, + 1 r(7-6i)(?i + i)i+» 



(29) 



so the initial slope is still (7 — 1)//?^ but there is a higher 
order term, which means that the plot of the ratio against 
1/(71 + 1) becomes curved. This behavior is simply a 
reflection of the true temperature dependence of xi^)- 

The second more annoying compfication consists of se- 
ries of "parasitic" terms with alternating signs which 
arise from the presence of anti- ferromagnetic poles 
Although when summed to inflnite n they give a zero or 
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negligible contribution to the true susceptibility xi^): in 
the ISG case they can lead to dramatic oscillations in the 
ratios a„+i/a„, Fig. [Ml and Fig. [511 The series can never- 
theless be analyzed, at least in dimension 4d and above, 
using the Fade approximant technique, accompanied by 
methods known as Ml and M2 0. 

An unorthodox but transparent variant on the graph- 
ical method is the following. Suppose the initial series 
for xiP) s-t some temperature /3 < /3c, terminating with 
term n, is written 

x) = 1 + aix + a2x'^ + a^x^ + ■ • ■ + a„x". (30) 

where x = . Then terms can be regrouped 

X*(n + l,x) = l+ {— + — \x 



4.T 2 ■ 

Ax 2 A 



032; \ 2 



X 



fan On+l CLn+2X\ ^+1 
^4^ 2 4 / ^ 

--l + alx^ h a;+ix"+^ 

= 1 + aix + • ■ • + a„x" 

4 4 



(31) 



The two sums x(n,x) and -I- l,a;) are identical up 
to the nth term. The series do not terminate in exactly 
the same place, but this is irrelevant as the aim is to 
extrapolate to infinite n so as to obtain the true total 
x(x). The essential point is that the sets of ratios of 
the successive coefficients R*{n + 1) = aj^+i/aj!! in the 
regrouped x* series now evolve smoothly with l/(n + 1), 
as can be seen in Figs. [26] and [27l and can be readily 
extrapolated to infinite n. This was certainly not the case 
for the raw series. The sum of the extrapolated regrouped 
terms can be considered a very good estimate of the true 
x{x). It turns out that for the 5d bimodal and Gaussian 
ISGs the set of regrouped term ratios are very insensitive 
to the choice of the trial and so the fit parameters for 
a representative provide good estimates for the true 
critical physical parameters. A simple fit R* {1 / {n+1)) = 
A + B{l/{n + 1) + C{l/{n + 1)^) provides estimates of 
the intercept A ~ l//3ci the initial slope B = (7 — 1)//?^ 
and the strength of the Wegner correction to scaling C = 
—a^9T{'y)/(T{'j — 9)B^) (For the fit we have assumed for 



convenience 
9 ). From Fig. 



^ 1 iTi but other values can be chosen for 



/3c = 0.3905, 7 = 1.85 and ~ 2.0 for 



the bimodal 5d ISG, and from Fig. ^ Pc = 0.420,7 = 
1.68 and ^ for the Gaussian 5d ISG. The values 
of the critical temperatures and the critical exponents 7 
are in excellent agreement with but appear to be more 
accurate than the central values from the much more 
sophisticated analysis of Daboul ei a/ Q- In addition, 
this method provides an estimate of the strength of the 



11.0 
10.5 
10.0 
9.5- 

_ 9.0 
c 

S. 8.5 



8.0 
7.5 
7.0 
6.5 
6.0 




0.4 0.6 
1/{n+1) 



FIG. 24: (Color online) The raw coefficient ratios a{n + 
l)/a{n) in the HTSE susceptibility series for the bimodal ISG 
in dimension 5 {a(n) data from [31). 
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FIG. 25: (Color onIine)The raw coefficient ratios a{n+l) / a{n) 
in the HTSE susceptibility series for the Gaussian ISG in 
dimension 5 (a(n) data from 0]). 



Wegner correction term, which was not explicitly cited 
as a result of the analysis in Ref. ■ 

It is important to underline that this technique is an 
analysis of the exact HTSE coefficients and so is entirely 
independent of the simulation data. Nevertheless the 
method again leads to a value for the bimodal critical 
exponent 7 ~ 1.85 which is quite different from the Gaus- 
sian critical exponent 7 ^ 1.65, in full agreement with the 
conclusions drawn from the simulation data. 
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FIG. 26: (Color online) The regrouped coefficient ratios 
a*(n + l)/a*{n) in the HTSE susceptibility series for the bi- 
modal ISG in dimension 5 for = 9. 
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FIG. 27: (Color online) The regrouped coefficient ratios 
a*(n+l)/tt*(n) in the HTSE susceptibility series for the Gaus- 
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